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We investigate low-energy properties of two-dimensional quantum spin systems with the ladder 
and plaquette structures, which are described by a generalized antiferromagnetic Heisenberg 
model with both of the bond and spin alternations. By exploiting a non-linear a model technique 
and a modified spin wave approach, we evaluate the spin gap and the spontaneous magnetization 
to discuss the quantum phase transition between the ordered and disordered states. We argue 
how the spin-gapped phase is driven to the antiferromagnetic phase in the phase diagram. 

KEYWORDS: 2D spin system, spin gap, non-linear a model, modified spin wave 



• §1. Introduction 



Recent intensive studies on high-T c superconductors 
and related materials have attracted renewed interest 
in the quantum phase transitions in low-dimensional 
spin systems. In particular, a variety of new materials 
have been found experimentally: For instance, CaV/iOg 
Jj],]^,||, which may be described by the two-dimensional 
(2D) Heisenberg model on a depleted square lattice, real- 
izes a quantum spin liquid even in two dimension. Also, 
SrCu2C>3 has been intensively studied as a proto- 

typical example of ladder compounds, ||,|7| for which the 
introduction of nonmagnetic impurities induces a novel 
phase transition to the magnetic state. For these ma- 
terials, the plaquette or ladder structure is essential to 
stabilize the non-magnetic phase with spin gap. The spin 
chains with an alternating array of different spins (mixed 
spin chains or alternating spin chains) also give another 
new paradigm, || [Tol [ll], |l2) for which the topology 
of spin arrangement plays an essential role to determine 
the low-energy properties. Quantum phase transitions 
in the above various low-dimensional spin systems have 
provided an interesting subject in quantum spin systems. 

In our previous paper, [jl3| by employing non-linear a 
model (NLcrM) techniques, we investigated the compe- 
tition between the gapful and gapless states in a gen- 
eralized spin ladder which possesses both of the spin 
and bond alternations. We systematically investigated 
not only ladder systems but also mixed spin chains and 
plaquette-type spin chains, and clarified how the com- 
peting interactions in addition to topological properties 
control the quantum phase transitions. In real mate- 
rials, however, the coupling between ladders or chains 
neglected in our previous paper becomes certainly im- 
portant, and may drive the system to the ordered state. 

Motivated by the above points, in this paper we study 
the quantum phase transitions between the ordered and 
disordered states in 2D spin systems, by extending the 
previous NLcrM analysis of the spin chain and ladder 
systems. In particular, we deal with the 2D antiferro- 
magnetic Heisenberg model with the ladder and plaque- 



tte structures, which also possesses the spatial variation 
of spins (mixed spins). By calculating the spin gap and 
the spontaneous magnetization by means of NLcrM and 
modified spin wave (MSW) approaches, |l4|, |l5| we dis- 
cuss how the disordered state with spin gap competes 
against the magnetically ordered state. 

The paper is organized as follows. In §2 we introduce 
a generalized 2D spin model for which both of the bond 
and spin alternations are included. In §3, by extending 
the method of Senechal, we calculate the spin gap 
at finite temperatures by means of NLcrM with saddle- 
point approximation, and then argue to what extent this 
approach can describe the quantum phase transitions for 
2D systems. We also point out the shortcoming in this 
approach for spin-1/2 systems. In § 4, we further study 
the model by using the MSW approach, which qualita- 
tively improves the results of NLcM in the spin-1/2 case. 
Brief summary is given in §5. 

§2. Model Hamiltonian 

We introduce a 2D antiferromagnetic spin system with 
both of the bond and spin alternations, which is de- 
scribed by the following Hamiltonian on a square lat- 
tice, 

H = J x [1 + (-1)*7«] ■ S i+1 j 

+ JvY / [l + (-l)hy]Si, j -S i , j+1 , (2.1) 

where J x (J y ) and j x (pj y ) are the exchange coupling and 
the bond-alternation parameter along the x (y) direction, 
respectively. Here Si j is the spin operator at the (i, j')-th 
site in the x — y plane. Since we deal with the antiferro- 
magnetic case with J x ,J y >0 and —1 < "fx-.Jy < 1, the 
ground state of the Hamiltonian (2.1) becomes the Neel 
state in the classical limit. What is distinct from an ordi- 
nary 2D square lattice model is that the present system 
includes not only the bond alternation but also the spin 
alternation. For simplicity, we here consider a specific 
arrangement of two types of spins shown schematically 
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Fig. 1. A generalized 2D spin system with both of the bond and 
spin alternations. 



in Fig. 1, where black circles and white circles denote 
the spins si and S2- The generalized Hamiltonian (2.1) 
allows us to study a variety of interesting spin systems 
which possess a ladder structure, a plaquette structure, 
and a mixture of different spins. This is a natural exten- 
sion of our previous analysis of the mixed-spin chains and 
ladders |l3| to more realistic cases including the effects 
of two dimensionality. In the next section we shall first 
study the model by employing the NLcrM techniques. 

§3. Non-Linear Sigma Model Analysis 

In this section we discuss the competition between the 
spin-gapped phase and the magnetically ordered phase 
by taking the non-linear a model (NLcrM) as a low en- 
ergy effective theory. Let us begin by briefly summariz- 
ing how NLcrM techniques are applied to our model in 
the continuum limit. Using the coherent-state path inte- 
gral formalism, |[7| the partition function in the present 
system is given by 



Z = J T)f2ijCxp — i SjjLu [f2jj] — J drH (t) 



(3.1) 



where 



w[f2ij]= / &TLp i3 , (1 - cosfly), (3.2) 
Jo 

H(t) = •/,. ^ S i,j S i • I ,' t 1 + ( _1 ) l 7x] ■ 

(3.3) 



In the above expressions we have introduced the 
unit vector fl by S = sf2, where fl = 
(sin 8 cos ip, sin sin ip, cos 8) . The term u> [fl] is the Berry 
phase. In the present approach we assume that the 
usual spin-wave analysis can correctly describe low- 
energy modes at wave vectors near (0, 0), (0, tt), (tt, 0), 
and (tt,tt). The last one corresponds to the ordering 
wave vector. Then the unit vector flij can be written 
in terms of three fluctuation fields and one staggered field 
as follows [fill 



«2i+a,2,-+/3 = (-l) a+f3 n(r') 



(3.4) 



n (r0 = 4> (r) + (-l) ma+nP tl (r) , (3.5) 

run 

where r' = r + am + (3n, m(n) and ot(/3) runs from to 
1, and a is the lattice constant. The primed sum means 
that the term m = n = is omitted. Therefore fluctua- 
tions which we take here are Iqi, ho, and hi for each site. 
The point we wish to notice is that the fluctuation I is 
further distinguished by the index k{— 1, 2), correspond- 
ing to spins si and S2, respectively. The staggered field 
4> is a slowly varying field on the scale of the lattice spac- 
ing. Substituting eq.(3.4) into eq.(3.3) and making the 
expansion up to quadratic order in I, <fi' , and </>, we arrive 
at the following Hamiltonian in the continuum limit 



H(t) = J dr(h x + h y ), 



(3.6) 



(4 + 4)(i + lx )(d x <t>y 



+2(l + 7 x)£W>- s&o 5 



h = ^ 



2(i-7,)E( s 1 4 1)2 + A ( ? )2 



(3.7) 



(1) , l(2) 
l 01 



(3.8) 



+ 4 SlS2 (l+ 7y ) [(d V <t>f + dyj> ■ (/,„ 

+ S i S2 (i+7,)E(4 1) +4 ) " 2 

i 

The Berry phase in this system is given by 

^Berry = i SjjW [flij\ 
ij 

= i J dr ( SiZ " + S2 '") ' x ^ (3 - 9) 

Integrating eq. (3.1) over the fluctuation fields I, we thus 
end up with the NLcrM 



Z = / Dc6cxp <^ - 



1 

~9 



v 



+vR x {d x 4>) 2 + vR y (d y </>) 2 ] } . (3.10) 



Note that the Berry phase terms cancel out in the present 
system and do not show up in the effective action. This 
simple result allows us to further investigate the quantum 
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phase transitions by employing ordinary NLcrM tech- 
niques. Since general forms of parameters g,v,R x , and 
R y are rather complicated, we shall show their concrete 
formulae for each case discussed in the following. 

3.1 Gap equation 

To proceed the calculation based on the NLcrM, fur- 
ther approximations may be needed. Here we use 
the saddle point approximation to obtain the spin gap 
equation from eq.(3.10), by extending the treatment of 
Senechal. |l(| The validity of this approximation will be 
discussed later in this section. In what follows, we set 
the speed v — 1; it can be restored by dimensional anal- 
ysis at the end of the calculation. Rescaling the field 
cf> — > ,/g(t>, the Lagrangian (3.10) is now written as 

l \(d T (t>) 2 +R X (d X Cj,) 2 +Ry(dy(t>) 2 



L = - 



-a {<f> 2 - i/g) 



(3.11) 



where a is the Lagrange multiplier to enforce the con- 
straint 4> 2 — 1/g. Integrating eq.(3.11) over the field 4>, 
we obtain the effective potential to the field a 

dk 



(27T) 



3.2 Plaquette- structure system 

Let us start our discussions for the spin system with 
plaquette structure with the parameters, 7 = 7, = 
7j/,J = J x = Jy, for which spins are assumed to be 
uniform (s = S\ = S2). In this case, the parameters in 
NLcrM are given by 

. 2\ 



R.r 
Ru 



1 



1 



= 2aVTTJ/s 1 



(3.14) 



2J y sa^/l +j, 



where j is defined by j — J x / ' J y . The cut-off parameter, 
as already mentioned, is taken as aA = \j2ir /3. 
The behavior of the normalized spin gap S(= A/Aq) 



CO 




0.1 0.2 



0.3 



0.4 0.5 



x log 1 + 



R x fc2 



Ryky 



(3.12) 



Here we replace the field a by the uniform value, which 
corresponds to the square of the spin gap. To evaluate 
the spin gap, the magnon speed v is restored by dimen- 
sional analysis. Finding out the saddle point of the po- 
tential eq.(3.12), we obtain the spin gap equation 

1 f dk 1 



fJ 



{2-Kf 



R x k x 



Ry ky 



(3 



x coth -\/R x kl 



R y kl + {A/vf. (3.13) 



It is to be noticed that a momentum cut-off parameter A, 
which is proportional to the inverse lattice spacing a , 
should be introduced in the above expression. Since it 
is hard to determine the cut-off procedure only from the 
present theory, we use other means to fix it. We outline 
a procedure for the case of plaquette structure charac- 
terized by the parameters: 7 = j x = j v ,J= J x = J y . 
In the limit 7 — > 1 the system changes to the assembly 
of plaquette singlets. Then the spin gap of the system 
is exactly A = 2J at T = 0. Requiring that the gap 
obtained from eq.(3.13) with 7 = 1 should be equal to 
2 J, the cut-off parameter is chosen as aA = ^/27r/3. As- 
suming that this value can be used for other values of 
7, we discuss the low-energy properties of spin systems 
with plaquette structure. 

This completes our formulation of the model in terms 
of NLcrM approach. In the following we show the nu- 
merical results obtained for the 2D spin systems with 
plaquette and ladder structures, as well as the mixed 
spin system. 
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Fig. 2. Normalized spin gap <5(= A/Ao) for spin s = 1/2 (a) and 
s = 1 (b) plaquette structures calculated by the NLcrM. It is 
shown as a function t(= aT). Aq is the spin gap at T = 0. 



is shown as a function of the temperature i(= aT) for 
spin s = 1/2 (Figs. 2(a)) and s = 1 (Figs. 2(b)) sys- 
tems. Ao is the spin gap at T = 0. Physically, the spin 
gap should be regarded as the inverse of the correlation 
length. As the temperature is decreased, the spin gap 
(inverse correlation length) decreases, and reaches zero 
or a finite value according to the nature of ground state. 
For the s — 1 case, we can clearly observe the quantum 
phase transition at T — in Fig. 2(b) when we decrease 
the value of 7 from unity. At 7 = 1 , the system at T = 
is composed of isolated plaquettes and forms the pla- 
quette singlet ground state which possesses the spin gap 
in the excitation. In decreasing 7, the plaquette-singlet 
state with spin gap is driven to the magnetically ordered 
state, which is realized for 7 < j c . As is the case for 
ordinary 2D quantum systems with continuous symme- 
try, the ordered state is allowed only for T = 0, which 
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is correctly described in our treatment, as is seen from 
the figure. Although the NLcM approach has produced 
qualitatively correct results for the case with s > 1/2, we 
can not find the quantum phase transition to the ordered 
state for s — 1/2, as seen from Fig. 2 (a). This is contra- 
dicted to the well-known fact that the ground state of the 
s = 1/2 isotropic square lattice (7 = 0) should be an an- 
tiferromagnetically ordered state. This implies that our 
approach may not sufficiently take into account antifer- 
romagnetic fluctuations, and this shortcoming may come 
mainly from the saddle point approximation we used. In 
order to improve our results for the s — 1/2 case, we fur- 
ther study the system by employing modified spin wave 
approach in the next section. 

3.3 Ladder- structure system 

Next we investigate the 2D system with ladder struc- 
ture which is realized by choosing j x = and 7 = j y . 
This is the system composed of periodic array of ladders 
along y-direction. p8| By appropriately tuning the in- 
teraction strength, we can naturally interpolate from the 
isolated ladders to the 2D system with ladder structure. 
The corresponding coupl ings for NLcrM is given by the 
same expression as (3.14). The cut-off parameter in the 
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Fig. 3. Normalized spin gap 5(= A/Ao) for the spin s = 1/2 (a) 
and s = 1 (b) ladder systems calculated by the NLerM (j = 1/2) . 
It is shown as a function t(= aT). 



ladder system is taken as aA — 0.99, which is determined 
according to the numerical results jl9],|2(]] (A/ J = 0.504) 
for the uniform isotropic ladder (j — 1/2). The results 
for the spin gap are shown in Figs. 3 for the cases of 
s = 1/2 and 1. Note that the system with the param- 
eter 7 = 0(1) corresponds to the isotropic square lat- 
tice (the isolated ladder). For the isolated ladder, the 



system should have the spin gap, which is consistent 
with our results. By decreasing 7 from unity, we ob- 
serve the quantum phase transition at T = from the 
spin-gapped phase to the magnetic phase for s = 1, but 
not for s = 1/2. As is the case for the model with pla- 
quette structure, we thus describe the quantum phase 
transition properly for s > 1/2. The problem we have 
encountered again for s — 1/2 has the same origin men- 
tioned before for the plaquette case. We shall come to 
this problem in the next section. 

3.4 Mixed spin system 

So far, we have considered the uniform-spin systems. 
We note that mixed spin chains with periodic array 
of several kind of different spins have also attracted 
much attention both experimentally and theoretically. 
|^,||;|l^,|ll[] For example, mixed spin chains with either 
the ferrimagnetic ground state or the singlet ground state 
have been intensively studied. |||l(|[Tl|] In the previous 
paper, |l3| ] we have studied the mixed spin chains by 
NLerM approach. To make the theoretical model more 
realistic, it may be natural to take into account the cou- 
pling among mixed chains. We address this problem in 
this section. As a simple example of mixed spin systems, 
we here consider the spin system composed of two dif- 
ferent spins Si = 1 and S2 = 3/2 (and also s\ = \/2 and 
S2 = 1) for the 2D system shown in Fig. 1. This type 
of system enables us to bridge 2D quantum systems and 
mixed spin chains mentioned above. To make our dis- 
cussions simpler, we deal with the system without bond 
alternation. Then parameters in NLerM for the mixed 
spin system are given by 



Rx 



j (1 + a) 2 (l + a 2 ) 



Ry = 1, 
1 



S2 

4JyS2a 
1 + a 



-aK, 



K, 



(3.15) 
(3.16) 
(3.17) 

(3.18) 



where K = y/(2j + 1)[(1 + a) 2 + 4ja]/[(l + a) 2 + 8ja] 
and a = S2/S1. The cut-off parameter is taken as aA 
1.11 according to the numerical results in ref pi]] : A/J — 
0.41 (T = 0, si — S2 — l,j = 0). The normalized spin 
gap calculated for si = 1 and S2 = 3/2 is shown in 
Fig. 4. For j = 0, the system is composed of isolated 
mixed spin chains, which is known to be massive, being 
consistent with the present results. |l3| When the value 
of j is increased, the system gradually changes to the 2D 
system. We then find a plausible result that there exists 
the critical value of j at which the system is changed to 
the ordered state. For the s± — 1/2 and S2 = 1 case , 
however, we have checked that in the whole range of j 
the system has the spin gap, and the long range order 
does not show up. Since we naturally expect that for 
moderately large j the antiferromagnetic ground state 
may be stabilized, we should improve our treatment for 
this case, which is reexamined in the next section by 
MSW approach. 
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Fig. 4. Normalized spin gap <5(= A/Aq) for the mixed spin sys- 
tem, for si = 1, S2 = 3/2, calculated by the NLtrM. It is shown 
as a function j(= Jx/Jy)- 



§4. Modified Spin- Wave Approach 

We have so far studied our 2D spin model by the 
NLcM with saddle point approximation, by employing 
the method of Senechal Jl6| , which is a natural extension 
of our previous study on the spin chain and spin ladder 
systems. 1 13 Although qualitative features for the quan- 
tum phase transitions have been described for the sys- 
tems with larger spins, we have encountered the problem 
for the s = 1/2 cases as well as the mixed spin case with 
s = 1/2,1, for which our NLerM analysis may not de- 
scribe the correct phase transition. We think that this 
shortcoming comes mainly ^from our saddle point ap- 
proximation, which may not properly incorporate mag- 
netic correlations for s = 1/2 case. To improve this 
point, there may be several possible ways. The one is to 
apply renormalization-group analysis to the NLerM in- 
stead of saddle point approximation, which was success- 
fully applied to 2D s = 1/2 systems. [^2[ Alternatively, 
there exists a different but similar approach based on 
the modified spin wave (MSW) approximation. p^,p^|] 
This method is known to provide rather satisfactory re- 
sults for the isotropic s — 1/2 spin systems, which are 
comparable to the renormalization group method men- 
tioned above. For instance, both of the above methods 
have produced the correct temperature dependence of 
the correlation length for the 2D s = 1/2 spin systems 
on square lattice. Therefore, the MSW approach is ex- 
pected to provide a sensible rout to improve our results 
for small s cases in the previous section. In what follows, 
we treat the present model within the MSW approxima- 
tion. We show the formulation in the case of uniform 
spins, for simplicity, since it is straightforward to extend 
it to the system composed of two distinct spins. 

By the Holstein-Primakoff transformation, the spin 
operators are expressed in terms of boson creation and 



annihilation operators as follows 



5+ = y 2s — ata r a r 
SI = blbr - s 



S r 

S r 
s r 



= bt\2s- btb r 

— C^* (s/j* & 



r — 



(2i, 2j) G A, 
(2z + l,2j) G B, 
(2i,2j + l) G C, 



S+ = V2s - drdr d r 



t. 



r = (2i+ l,2j + l)eD. 



(4.1) 

Boson operators a r , b r , cy, d r belong to the A, B, C, D 
sublattices, respectively. Since the unit cell in our sys- 
tem includes four sites, we have introduced four kinds 
of Bose operators. In case that we consider the system 
with two kinds of spins, eight kinds of bosons should be 
introduced. The Fourier transformations are defined as 




for a — a, d(/3 — b,c), where N is the total number of lat- 
tice, and the summation over momenta k is restricted to 
the quarter of the original Brillouin zone. Then Hamil- 
tonian becomes 

H = 2 J v s { (1 + 3) + b t bk + c t Ck + d t dk ) 

k 

+jK x Ua k b k + 44) e' 01 ' + H 

[( 



[ a k c k + 44) e " 



H.c. 



} 



(4.3) 



where 



k x = \Jl - (1 - 7|) sin 2 k x , 
4> x = arg [cos k x + i-y x sin k x ] 



K y = yjl - (1 - 7y)sin 2 k y , 

tfiy = arg [cos k y + vy v sin k y ] . 

The diagonalization of the Hamiltonian (4.3) yields the 
conventional spin wave excitations. The staggered mag- 
netization per site is expressed as 



1 

N 



EOZ \ ^ QZ \ ^ QZ 

D r 2-^1 r 2-^1 r 



.r£A 



r£B 



r£C 



r£D 



jj ]T [ajofc + b\b k + c\c k + dld k ] . (4.4) 



Here it should be noticed that the sublattice symmetry 
Jig ] in the ground state is broken. In order to extend our 
study into the region of the disordered phase, we restore 
this symmetry by enforcing the following constraint on 
eq-(4.3), 



(4.5) 



This means that the staggered magnetization (4.4) 
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should be zero in the disordered phase. Since the con- 
straint may give rise to the spin gap in a certain param- 
eter region, we can discuss the phase transition ^from 
the ordered phase to the disordered phase with spin gap. 
Using the Lagrange multiplier A, the Hamiltonian with 
the constraint (4.5) is cast to 



H = H -A 



E [ a l ak 



btbk 



d\d k 



• (4-6) 



By the Bogoliubov transformations 



i 

V d k ) 



1 

71 



v k + 



Vu 



where 



sinh^e 1 ^ 



cosh^j. 



4 2)t 
V /f )f J 



sinhfl^e"^ 
— cosh 9^ 



(4.7) 



-sinh^ 



cosh (9^^ 



cosh (9^ e-^ 



coshfl., = —=a / — r 




M 1 2(J X + J v )s' 
the Hamiltonian (4.6) is diagonalized as 



(4.8) 
(4.9) 

(4.10) 

(4.11) 
(4.12) 



H = 



EEK^^+^^Vf]. (4.13) 



k i=l 



where the energy dispersion E k is defined by 2(J X 



J y )sfjLe^. The spin gap A opens at fc 
estimated as 



A = J\ 2 - 4(J X + J y )sX. 



(0, 0), and is 



(4.14) 



We can easily confirm that in the case of A = this the- 
ory is identical with the conventional spin-wave theory 
and that the lowest mode of the system has the linear 
dispersion. 

In order to formulate thermodynamics, we introduce 
the probability, P^.(nf) {P^ k {n~)) 1 for the state with the 
wave vector fc being occupied by nf(n~) bosons with 
type aW(/? (l) ). Then the free energy has the form 



F = 



k,i,n 



{EinP+{n)+E^nP- k {n) 



[^tk P Un)+^Pr k {n)]}-TS, (4.15) 
where fj,f k is the generalized chemical potential, and S is 



the entropy. By evaluating Pju(ri) which minimizes the 



free energy F, we end up with an ordinary formula, 

P± = (1 - exp [-E±/T]) exp [-nE±/T] . (4.16) 



Taking the N 
form 

I 7T Z 



co, the resulting constraint has the 



dfc„ 



1 * E k 1 u E k 
— coth-^ + — coth-^ 



(4.17) 



This completes the MSW formulation of our model. 
Solving the above equation to determine the chemical 
potential A, we the n eva luate the spin gap at finite tem- 
peratures from eq.( 4.14 ). 

In the following we discuss the results separately for 
three cases studied in the previous section. 

4-1 Plaquette- structure system 

Let us first show the results obtained for the system 
with plaquette structure. We here focus on the s = 1/2 
system for which our NLctM approach with saddle-point 
approximation has encountered the problem. The present 

1.4 
1.2 

1 

< °- 8 

0.6 

. 4 
0.2 



1=1 













0.2 



. 4 



. 6 



. 8 



Fig. 5. Spin gap A for the plaquette structure calculated within 
the MSW approximation. It is shown as a function of the 
temperature? 1 . 
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Fig. 6. Staggered magnetization m (the dashed line) and the spin 
gap A (the solid line) for the plaquette structure calculated 
within the MSW approximation. 



MSW approach provides us with not only the tempera- 
ture dependence of the spin gap (Fig. 5) but also the 
staggered magnetization at T = (Fig. 6). We find in 
Fig. 5 that, in contrast with the results of the NLerM, the 
solid line for 7 = reproduces the gapless behavior at 
T = 0. Moreover we can see in Fig. 6 that when 7 = 0, 
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the system has a finite magnetization (to = 0.303). The 
magnetization decreases with increasing 7 and finally 
vanishes at the critical value 7 C = 0.798. For 7 > j c 
the plaquette-singlet ground state with spin gap is stabi- 
lized, as mentioned in the previous section. These results 
qualitatively improve the results of NLerM with saddle- 
point approximation since the effect of antiferromagnetic 
correlation may be treated more properly in the MSW 
approximation. Though the MSW can qualitatively de- 
scribe the correct behavior of the quantum phase tran- 
sition, it may be plausible to evaluate various quantities 
more quantitatively To this end, we have also performed 
a complementary calculation based on the cluster ex- 
pansion starting from the isolated spin plaquette, up to 
fourth order, which has been developed by several groups 
recently. J23| By applying the Pade approximation to the 
above series expansion, we have confirmed that there in- 
deed exists the quantum phase transition and the critical 
value 7 C is around 0.3. |24| This largely improves the re- 
sult of the MSW quantitatively. The detail for the results 
of the cluster expansion will be reported elsewhere. 

4-2 Ladder- structure system 

We next consider the 2D s = 1/2 system with ladder 
structure discussed in the previous section. |l§| l In Fig. 7, 
the spin gap calculated by the MSW is shown as a func- 
tion of the temperature for several choices of the bond 
alternation parameter. It is seen from Fig. 7 that there 
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massive, and thus we encounter the second phase transi- 
tion. 
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Fig. 8. Phase diagram of the ladder system calculated by the 
MSW approximation (see the results of ref. [jl^]). The broken 
line is the phase boundary determined by the dimer expansion 
method up to fifth order combined with the Pade approximants. 
We have also shown the results of the Quantum Monte Carlo 
calculation |18[ by the dots. 



We should note here that the present ladder-structure 
system at T = was previously studied by Katoh and 
Imada by MSW theory, p) Our results at T = in- 
deed coincide with their results. Furthermore, they have 
performed the Monte Carlo calculations, and confirmed 
the quantum phase transition at the critical strength of 
the bond-alternation parameter, and estimated its value 
rather precisely, which improves the results of the MSW 
quantitatively. In this connection, we have also per- 
formed the calculation based on the dimer expansion ]23j ] 
up to fifth order by taking isolated dimers for rungs as 
a starting point. This approach may be complemen- 
tary to the Monte Carlo approach of Katoh and Imada. 
The phase boundary determined by the dimer expansion 
method with Pade approximation is plotted in Fig. 8 
with broken line for smaller j region, which shows fairly 
good agreement with Monte Carlo results. 



Fig. 7. Spin gap A in the s = 1/2 ladder structure for j = 2 
calculated by the MSW approximation. 

is a critical strength of the bond-alternation parameter 
7 C at which the correlation length (inverse of the spin 
gap) becomes large with decreasing temperature, and di- 
verges at zero temperature. For 7 < 7 C , the spontaneous 
antiferromagnetic order emerges, but only at zero tem- 
perature. In Fig. 8, we show the phase diagram in terms 
of the bond alternation parameter 7 and the anisotropic 
parameter j = J x /J y . If we fix 7 to be an appropriate 
value, the quantum phase transitions may be observed 
twice with the increase of j. If j is very small, the sys- 
tem is approximated by weekly coupled spin chains with 
dimcrization. This phase is massive because of the dimer 
gap for the spin chain. When j is increased, the system 
is driven to the 2D antiferromagnetically ordered state. 
If j is further increased, the system may be described by 
a weakly coupled ladder system which should be again 



4-3 Mixed spin system 

As mentioned in the previous section, the mixed spin 
chains have attracted much current interest. Although 
we could describe the quantum phase transition for the 
s = 1,3/2 mixed system qualitatively by NLcrM ap- 
proach, we failed to do for the s = 1/2,1 mixed case, 
which may be most interesting from the viewpoint of 
experiments. We address the latter case in this subsec- 
tion. In Fig. 9, we show the results for the spin gap and 
spontaneous magnetization calculated by MSW method 
for mixed spin system with s = 1/2, 1. For j = 0, the 
system is reduced to the mixed spin chain with the pe- 
riodic array of spins as • • • ll^ll^j • • •, so that it is 
in the massive phase because of the topological nature 
of the system. jl3| This is consistent with the results 
of NLerM. In contrast to NLerM approach, however, the 
spin gap rapidly decreases in the vicinity of j — 0, and 
the systems is driven to the ordered phase with a finite 
staggered magnetization at j c = 0.009. 
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Fig. 9. The staggered magnetization m (the dashed line) and the 
spin gap A (the solid line) for the s = 1/2, 1 mixed spin system 
calculated by the MSW approximation. The calculated quanti- 
ties are shown as a function j(= J x /Jy)- 



Up to now, the mixed spin chains found experimentally 
are known to exhibit the ferrimagnetic ground state. ||] 
It may be expected that the mixed spin chains with sin- 
glet ground state, as discussed here, may be also fab- 
ricated experimentally in the future. Then such mixed 
spin chains should provide an interesting paradigm of 
quantum spin systems, for which the topological nature 
of spatial variation of spins control the low-energy prop- 
erties of the systems. In particular, it may be interesting 
to observe how such massive systems are driven to the 
magnetically ordered phase in the presence of the inter- 
chain coupling, impurities, etc. 

§5. Summary 

By extending our previous studies on the spin chain 
and ladders, we have systematically investigated quan- 
tum phase transitions in 2D spin systems with plaque- 
tte and ladder structures. In these systems, the lat- 
tice structure as well as the competing interactions play 
a crucial role to determine low-energy properties. By 
computing the spin gap and the spontaneous magnetiza- 
tion, by NLerM and MSW approaches, we have confirmed 
that both of the above approaches describe the quantum 
phase transition qualitatively well, so far as the system 
with larger spins is concerned. However, it has turned 
out that for s = 1/2 systems (as well as s — 1/2, 1 mixed 
spin systems), the NLerM with saddle point approxima- 
tion cannot properly describe the phase transition. So, 
this approach employed recently by Senechal may not be 
appropriate to discuss the phase transition even at qual- 
itative level for s = 1/2 spin models. On the other hand, 
the MSW treatment, which may be provide the results 
comparable to the renormalization group treatment of 
NLerM, has been shown to plausibly describe the phase 
transition for s — 1/2. 

In order to give more quantitative discussions for the 
quantum phase transitions, we have also shown some 
preliminary calculations based on the series expansion 
methods such as the dimer expansion, the plaquette ex- 
pansion, etc, which indeed improve the results of NLerM 
and MSW. To derive more precise information on these 
quantum spin systems, it is desirable to extend the clus- 



ter expansion to higher orders systematically, which is 
now under consideration. 
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